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We regard the Wheeler-De Witt equation as a Sturm-Liouville problem with the cosmological 
constant considered as the associated eigenvalue. The used method to study such a problem is a 
variational approach with Gaussian trial wave functionals. We approximate the equation to one loop 
in a Schwarzschild background. A zeta function regularization is involved to handle with divergences. 
The regularization is closely related to the subtraction procedure appearing in the computation of 
Casimir energy in a curved background. A renormalization procedure is introduced to remove the 
infinities together with a renormalization group equation. 
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I. INTRODUCTION 

One of the most fascinating and unsolved problems of 
the theoretical physics of our century is the cosmological 
constant. Einstein introduced his cosmological constant 
A c in an attempt to generalize his original field equations. 
The modified field equations are 

R^v - ■zdnfR + ^c9^u = SttGT^, (1) 

where A c is the cosmological constant, G is the gravita- 
tional constant and 2),„ is the energy-momentum tensor. 
By redefining 



rptot _ rp 

1 fj,u — 



8ttG 



(2) 



one can regain the original form of the field equations 
R»» ~ \g^R = 8ttGT*°* = SttG (v + l£) , (3) 

at the prize of introducing a vacuum energy density and 
vacuum stress-energy tensor 



PA 



8ttG' 



rpA _ 



Alternatively, Eq.QJ can be cast into the form, 



R^v - -^g^uR + Aeffguv — 0, 



(4) 



(•5) 



where we have included the contribution of the vacuum 
energy density in the form T^ v = — (p) g^. In this case 
A c can be considered as the bare cosmological constant 



A e/ / = 8irGp ef f = A c + 8ttG (p) 



(0) 



Experimentally, we know that the effective energy den- 
sity of the universe p e ff is of the order 10 'GeV . A 
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crude estimate of the Zero Point Energy (ZPE) of some 
field of mass m with a cutoff at the Planck scale gives 



E 



ZPE 



1 

2j» 



d 3 k 
(2^ 



Vk 2 



16tt 2 



10 71 GeV 4 



(7) 

This gives a difference of about 118 orders [ij. The ap- 
proach to quantization of general relativity based on the 
following set of equations 



2k 

71' 



^(R-2A C ) 
Zk 



and 



-2V<7r«*k«]=0, 



* \9ij\ = (8) 



(9) 



where R is the three-scalar curvature, A c is the bare cos- 
mological constant and k = 8irG, is known as Wheeler- 
De Witt equation (WDW)0. Eqs.© and © describe 
the wave function of the universe. The WDW equation 
represents invariance under time reparametrization in an 
operatorial form, while Eq.@ represents invariance un- 
der diffcomorphism. Gijki is the supermetric defined as 

Gijki = ^(gikgji + giigjk - g^gu)- (10) 

Note that the WDW equation can be cast into the form 



2k 

Vg 



Gi-jkin^n 



2k 



R 



Vg 



A c * [gi. 



(ii) 

which formally looks like an eigenvalue equation. In this 
paper, we would like to use the Wheeler-De Witt (WDW) 
equation to estimate (p) . In particular, we will com- 
pute the gravitons ZPE propagating on the Schwarzschild 
background. This choice is dictated by considering that 
the Schwarzschild solution represents the only non-trivial 
static spherical symmetric solution of the Vacuum Ein- 
stein equations. Therefore, in this context the ZPE can 
be attributed only to quantum fluctuations. The used 
method will be a variational approach applied on gaus- 
sian wave functional. The rest of the paper is structured 
as follows, in section [H] we show how to apply the vari- 
ational approach to the Wheeler-De Witt equation and 
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we give some of the basic rules to solve such an equation 
approximated to second order in metric perturbation, in 
section ITTT1 we analyze the spin-2 operator or the opera- 



tor acting on transverse traceless tensors, in section lTlI Al 
we analyze the spin-0 operator, in section llVl we use the 
zeta function to regularize the divergences coming from 
the evaluation of the ZPE for TT tensors and we write 
the renormalization group equation, in section IIV Al we 
use the same procedure of section llVl but for the trace 
part. We summarize and conclude in section IVI 



II. THE WHEELER-DE WITT EQUATION AND 
THE COSMOLOGICAL CONSTANT 



where h is the trace of hij and is the three dimen- 
sional scalar curvature. To explicitly make calculations, 
we need an orthogonal decomposition for both mj and hij 
to disentangle gauge modes from physical deformations. 
We define the inner product 



(h, k) := / ^G^ kl h l3 (x) k kl (x) d 3 x, (18) 

JY, 

by means of the inverse WDW metric Gijki, to have a 
metric on the space of deformations, i.e. a quadratic 
form on the tangent space at hij , with 



The WDW equation © , written as an eigenvalue equa- 
tion, can be cast into the form 



A E * \gij] = A'* [ gi j} , 



(12) 



Qijkl = ( g ik gJ l + g il g jk _ 2g l ig kl ). 



(19) 



The inverse metric is defined on cotangent space and it 
assumes the form 



where 



A' = —— 



(13) 



(P, l) ■= / V9GijMP ij (x) q kl (x) d 3 x, (20) 



so that 



We, now multiply Ea. H12|) by ^* [g^] and we functionally 
integrate over the three spatial metric <7ij, then after an 
integration over the hypersurface S, one can formally re- 
write the WDW equation as 



kd 3 xA 1 



V 



mm 



V SV[gij}**[gij]*[ gi j] 

(14) 

The formal eigenvalue equation is a simple manipulation 
of Eq.JSJ. However, we gain more information if we con- 
sider a separation of the spatial part of the metric into a 
background term, g^ , and a perturbation, hij , 



gij gij ~t~ hi 



Thus ea. lfHjl becomes 



A 



(o) 



A 



(i) 



<*!*> 



where A^ represents the i th order of perturbation in h^ . 
By observing that the kinetic part of As is quadratic in 
the momenta, we only need to expand the three-scalar 
curvature J d 3 x^/gR^ up to quadratic order and we get 



(15) 



= A'* [g^ , 
(16) 



/ d 3 xy/g 



-hAh + -h u Ah H - -h ij ViVM 
4 4 2 3 



+ \hV l V l h u - ~h ij Ri a h1 + \hR l0 h 13 + ~h (i? (0) ) h 

(17) 



£iijnrn q 



(21) 



Note that in this scheme the "inverse metric" is actu- 
ally the WDW metric defined on phase space. The de- 
sired decomposition on the tangent space of 3-metric 
=deVbrmationsl3lllEli is: 



(22) 



Mj = 2 h 9H + ( L 0ij + Kj 



where the operator L maps £j into symmetric tracefree 
tensors 



(LOi; = + V,C - - 9i j (V • • 



(23) 



Thus the inner product between three-geometries be- 
comes 

(h, h) ■= J JgG ijkl hij {x) h k i (x) d 3 x = 



V~9 



- h 2 + (L0 ij (L0ij + h ij± hij 



(24) 



With the orthogonal decomposition in hand we can define 
the trial wave functional as 

* [hij (1?)] = AAl> [h^ (^)] * [hi (1?)] * [ h %™ ("£)] , 

(25) 
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where 



[^(^)]=exp{-i<Wf- 1 < B } 

=exp{-i<(^)^- 1 (^))! i J ■ (26) 

* ( - )] = ex P {-i <^-^):;; ce } 

The symbol "_L" denotes the transverse-traceless tensor 
(TT) (spin 2) of the perturbation, while the symbol "||" 
denotes the longitudinal part (spin 1) of the perturba- 
tion. Finally, the symbol "trace" denotes the scalar part 
of the perturbation. J\f is a normalization factor, (•, ■) x 
denotes space integration and K^ 1 is the inverse "prop- 
agator" . We will fix our attention to the TT tensor sec- 
tor of the perturbation representing the graviton and the 
scalar sector. Therefore, representation <|25|) reduces to 



-j\f 



exp {-1 (hK-'h)^ exp j-± (hK^h) 



Trace 
x,y 



(27) 

Actually there is no reason to neglect longitudinal per- 
turbations. However, following the analysis of Mazur 
and Mottola of Ref.0 on the perturbation decomposi- 
tion, we can discover that the relevant components can 
be restricted to the TT modes and to the trace modes. 
Moreover, for certain backgrounds TT tensors can be a 
source of instability as shown in Refs. [7j • Even the trace 
part can be regarded as a source of instability. Indeed 
this is usually termed conformal instability. The appear- 
ance of an instability on the TT modes is known as non 
conformal instability. This means that does not exist a 
gauge choice that can eliminate negative modes. To pro- 
ceed with Eq. H16[) . we need to know the action of some 
basic operators on \& [hij]. The action of the operator hij 
on |*) = \fr [hij] is realized by@ 



hij(x)\y) = h^ (~x) * [hij] . 



(28) 



The action of the operator tt^ on \*5f), in general, is 



be stationary against arbitrary variations of VP [hij]. Note 
that Eq.J2U can be considered as the variational analog 
of a Sturm-Liouville problem with the cosmological con- 
stant regarded as the associated eigenvalue. Therefore 
the solution of Ea. ljl4fl corresponds to the minimum of 



Eq.dSD. Thc form of (* 



A, 



\1/ ) can be computed with 



the help of the wave functional Ij27(l and with the help of 





(*N (*)!*) 



<*!*> 



(32) 



Kijki {~x, ~v) 



Since the wave functional l|27|) separates the degrees of 
freedom, we assume that 



A' = A 
then Ea. l|31|) becomes 

i (*|A 



A' 



v (*[*) 
WW) 



A 



= A" 



'_L 



(33) 



(34) 



Extracting the TT tensor contribution, we get 
1 
W 



dPxy/]jG l3kl (2k)K~ 1x {x,x 



ijkl 



(35) 



The propagator K (x,x) iakl can be represented as 



2A(r) 



where h['^ >± (~x) are the eigenfunctions of A2. t denotes 
a complete set of indices and A (r) are a set of varia- 
tional parameters to be determined by the minimiza- 
tion of Ea. (|35|l . The expectation value of A^ is easily 
obtained by inserting the form of the propagator into 
Eq.CH> 

A' x (A 4 ) 



TTij (x) I*) 



Shij (x) 



* [hij] , 



(29) 



while the inner product is defined by the functional inte- 
gration: 



(* 1 I * 2 ) 
We demand that 



[Vhij] *I [hij] * 2 [h k 



kl 



(30) 



1 



2 r 



[EE 

t i—1 



(2k) X z (t) + 



(2k) A, (r) 



(37) 



By minimizing with respect to the variational function 
Xi (t) , we obtain the total one loop energy density for 
TT tensors 



(38) 



V 



(31) 



where 

A tt (\ 1 )=^ ± (\ 1 ) = -A ± (\ 1 )/k. (39) 
The above expression makes sense only for u>f (r) > 0. 
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III. THE TRANSVERSE TRACELESS (TT) 
SPIN 2 OPERATOR FOR THE 
SCHWARZSCHILD METRIC AND THE W.K.B. 
APPROXIMATION 



From the transversality condition we obtain h\ (r) = 
h\ (r) . Then K (r) = L (r) . For a generic value of the an- 
gular momentum L, representation (|46[1 joined to Ea. 142(1 
lead to the following system of PDE's 



The Spin-two operator for the Schwarzschild metric is 
defined by 

(A 2 h TT Y. := - (A T h TT ) j . + 2 {Rh TT ) j . , (40) 

where the transverse-traceless (TT) tensor for the 
quantum fluctuation is obtained by the following 
decomposition 



(41) 



This implies that (fr T )^ Sj = 0. The transversality con- 
dition is applied on (h T y. and becomes Vj (h T Y- = 0. 
Thus 



(A T h 



TT\3 



-*s(^)i + Ul-^), (42, 



where As is the scalar curved Laplacian, whose form is 

\ - ( 1 2MG \ d2 + / 2?, -3MG\ d L 2 
\ r J dr 2 \ r 2 J dr r 2 

and R°j is the mixed Ricci tensor whose components are: 



R" 



2MG MG MGj 



(44) 



This implies that the scalar curvature is traceless. We are 
therefore led to study the following eigenvalue equation 



{A 2 h™)\=^h 



2 h i 

j 



(45) 



where lo 2 is the eigenvalue of the corresponding equa- 
tion. In doing so, we follow Regge and Wheeler in analyz- 
ing the equation as modes of definite frequency, angular 
momentum and parity 0. In particular, our choice for 
the three-dimensional gravitational perturbation is rep- 
resented by its even-parity form 



with 



diag [H (r) ,K(r),L (r)] Y lm (0, 



H(r) = h\ (r)-ift(r) 
K(r) = h 2 (r)-U(r) 
L{r)=h\ (r)-|A(r) 



(46) 



(47) 



-As 



(1 



2MG 



2MG 



)-*%£)H(r)=u* l H(r) 



Defining reduced fields 



H(r) = 



/iW. 



K{r) = 



h{r) 



(48) 



(49) 



and passing to the proper geodesic distance from the 
throat of the bridge 

dr 



dx = ± 



the system (I48|l becomes 



1 - 



2A/G 



(50) 



-& + Vi(rj\fi(x)=<4,f 1 (x) 



(51; 



-£ + V 2 (r) f 2 {x)=u J 2 J 2 {x) 



with 



(52) 



(53) 



Vi (r) = ii^l + Ui (r) 

V 2 (r) = + U 2 (r) ' 
where we have defined r = r (x) and 

f/i(r) = [^(l-^)-^] 

Note that 

U x (r) > when r > 
U x (r) < when 2MG < r < 

U 2 (r) > Vr e [2MG, +oo) 

The functions U\ (r) and U 2 (r) play the role of two r- 
dependent effective masses m 2 (r) and m\ (r), respec- 
tively. In order to use the WKB approximation, we define 
two r-dependent radial wave numbers k\ (r,l,u>i. n i) and 
k 2 (r,l,oj 2>n i) 



(54) 



k 2 (r,I, Wl , nI )=w?, r u- i ^-"»? M 
k 2 (r,l,tj 2 , nl )=ul nl -^-m 2 (r) 



(55) 



for r > 5^2. When 2MG < r < ^p, k 2 [r,l,Lu hnl ) 
becomes 

u2 r„ , . . x . .2 1(1 + 1) , _ 2 



+ (r) . (56) 
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A. The trace part contribution 

The trace part of the perturbation can be extracted 
from Ea. (|17fl to give 



d 3 xy/g 



-hAh + -hVi\/ih li 



--hR in h i] +~h(RW) h 



(57) 



(58) 



Some remark before proceeding is in order. First of all 
one has to note the minus sign in front of the laplacian 
acting on the scalar term: this is the origin of the con- 
formal instability. Secondly the use of gaussian wave 
functionals leads to a further simplification in Eci. (|58J) . 
because mixed terms of the type h (some operator) h TT 
disappear. Last for some backgrounds like the one we 
have considered here the three scalar curvature vanishes. 
Therefore the only piece of the quadratic order coming 
from the scalar curvature expansion is 



/ d 3 xy/g 



-\h/\h 
4 



(59) 



Moreover, if we follow the orthogonal decomposition of 
Eci. (|22[) . we can write 



1 



lonq . 



then the trace part of As becomes 
- 2 1 



2k 
~6 



drx— 
s V9 



2 K j/ X ^{-t Ah 



(60) 



(61) 



By repeating the scheme of calculation of section[n]to the 
trace part in Eq. I|61() , we get the scalar part contribution 
to the ZPE 



A" trace 
S 



\{2k)K~ 1 {x,x) 
6 



6 (2k) A (r) 



(64) 



By minimizing with respect to the variational function 
A (r), we obtain the total one loop energy density for the 
scalar component 



A' 



(A) 



(65) 



where A trace (A) = (A') (A) = -A*™ ce (A) /k. The 
above expression makes sense only for u> 2 (r) > 0. In 
the Schwarzschild background, the operator A can be 
identified with the operator As of Ea. l|43() . If we repeat 
the same steps leading to Ea. lfBT)) . we get 



d? 
dx 



—2+V{r) 



where we have defined 



h(r) 



f(x)=u 2 f(x) 



f(r) 



(66) 



(67) 



and used the proper distance from the throat defined by 
Ea.l |50|) . In order to use a W.K.B. approximation, we 
define a r-dependent radial wave number k(r,l, uj n i) 



k 2 (r,l,uj nl ) = uj 2 nl 



1(1 + 1) MG 



(68) 



IV. ONE LOOP ENERGY REGULARIZATION 
AND RENORMALIZATION 

In this section, we proceed to evaluate Ea. l|38(l . The 
method is equivalent to the scattering phase shift method 
and to the same method used to compute the entropy in 
the brick wall model. We begin by counting the number 
of modes with frequency less than uii, i = 1,2. This is 
given approximately by 



(2k 



(A) K (x, x) 



The propagator K (x, x) can be represented as 

(-x)h^ (~y) 



K 1 - 0?,7) :=J2- 



2A(r) 



(62) 



(63) 



where h^ T ' (~x) are the eigenfunctions of A. r denotes 
a complete set of indices and A (r) are a set of vari- 
ational parameters to be determined by the minimiza- 
tion of Eq. H62|) • The expectation value of A|T ace is easily 
obtained by inserting the form of the propagator into 
Eq.® 



A 



'trace 



(A,) 



9 (wi) 



Ui(!,Wi) (21 + 1). 



(69) 



where vi (l,uji), i = 1,2 is the number of nodes in the 
mode with (l,cJi), such that (r = r (x)) 



1 f+°° I 

v i{ l ,Ui) = —l dxJk 2 (r,l,Ui). 
2tt J^oo v 



(70) 



Here it is understood that the integration with respect 
to x and / is taken over those values which satisfy 
kf (r,l,Wi) > 0, i = 1,2. With the help of Eas. (l69l 70ll . 
we obtain the one loop total energy for TT tensors 



if 

8tt ^ 



+oo 



dx 



+oc dg (uji) 

Wi : dWi 

diO; 



(71) 
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By extracting the energy density contributing to the cos- 
mological constant, we get 



A TT = 



1 



16tt 2 



+ 00 



J?y/u? + m 2 (r)d(Ji, (72) 



where we have included an additional 47r coming from the 
angular integration. We use the zeta function regulariza- 
tion method to compute the energy densities p\ and pi. 
Note that this procedure is completely equivalent to the 
subtraction procedure of the Casimir energy computa- 
tion where the zero point energy (ZPE) in different back- 
grounds with the same asymptotic properties is involved. 
To this purpose, we introduce the additional mass pa- 
rameter /x in order to restore the correct dimension for 
the regularized quantities. Such an arbitrary mass scale 
emerges unavoidably in any regularization scheme. Then 
we have 



Pi (e) 
where 



1 



,2 E 



16tt 2 



-Hoc. 



dujj 



{(jf + rn 2 (r)) £ 



(73) 



Pi 



P2 



16 



h Jo °° w iV w i + TO ? (r)dui 



(74) 



16 



h Jo °° w 2 V^l + m l {r)dui 



The integration has to be meant in the range where usf - 
m 2 (r) > 1 . One gets 



Pi (e) 



256tt 2 



£ 



m| (r) 



2 In 2 - - 
2 



(75) 



-m% (r) 



In 



ml (r) 



21n2 - - 
2 



The quantity in Eq.JTUJ depends on the arbitrary mass 
scale /i. It is appropriate to use the rcnormalization group 
equation to eliminate such a dependence. To this aim, 
we impose that[lO| 



1 flAgrfr) d 



8ttG 



(80) 



Solving it we find that the renormalized constant Aq T 
should be treated as a running one in the sense that it 
varies provided that the scale /i is changing 

A™ ( M , r) = \r W r) + ^~ K M + ^ (0) In ^. 

1D7T /Xo 

(81) 

Substituting Eq.(5TJ into Ea.lfftty we find 



8ttG 



256tt 2 



ml W 



In 



mf (r) 
Mo 



-m£ (r) 



'"'IF)- 21 " 2 + 5 



- 21n2+ - 
(82) 



(83) 



In order to fix the dependence of A on r and M, we find 
the minimum of A$ T (/iq , r) . To this aim, last equation 
can be cast into the form 2 



i = 1,2. In order to renormalize the divergent ZPE, we 
observe that from Ea. H72f> . after reinserting the gravita- 
tional constant, we write 



ATT 
Jin 



(Po,r) 



l4 



8irG 



256tt 2 



x 2 (r) 



In 



\x(r)\ 



(84) 



A TT = -8^G ( Pl (e)+p 2 (e)). 



(76) 



To handle with the divergent energy density we extract 
the divergent part of A TT , in the limit e — > and we set 



A 



TT.div 



G 
32ire 



(mf (r) 



W) 



(77) 



+2/ 2 M 



In 



2/M 



(85) 



where a; (r) = ±to 2 (r) //x 2 , and y (r) = m 2 (r) //i 2 ,. Now 
we find the extrema of Ao (/xo; x (r) ,y (r)) in the range 
< r and we get 



Thus, the renormalization is performed via the absorp- 
tion of the divergent part into the re-definition of the 
bare classical constant A TT 



A 



TT 



Aq T - A 



TT,< 



(78) 



The remaining finite value for the cosmological constant 
reads 



x(r) = 
y (r) = 



which is never satisfied and 



x (r) = 4/e 
y(r) = 4/e 



(86) 



m\ (r) = 4xi 2 /e . . 
m 2 (r)=4^ 2 /e " ^ 



1 



AT = 

8ttG 256tt 2 



(r) 



In 



m 2 (r) 



21n2- - 
2 



1 Details of the calculation can be found in the Appendix. 



2 Recall Eqs. I55I 561 . showing a change of sign in m\ (r). Even if 
this is not the most appropriate notation to indicate a change 
of sign in a quantity looking like a "square effective mass" , this 
reveals useful in the zeta function regularization and in the serch 
for extrema. 
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which implies M = and r — \/3e/2/Xo- On the other 



hand, in the range IMG < r < 5A ^ G , we get again 



x (r) = 
«(r)=0 ' 



(88) 



Thus, the renormalization is performed via the absorp- 
tion of the divergent part into the re-definition of the 
bare classical constant A trace 



A' 



p^trace ^trace,div 



(97) 



which has no solution and 



which implies 



—m\ (r) — 
ml (r) = 4/ig/e 



(89) 



Therefore, the remaining finite value for the cosmological 
constant reads 

A trace 

8ttG 



M = 4^ f 3/ 3eG 
f=V6e/Vo 
Eg. 1)85(1 evaluated on the minimum, now becomes 



G 



2e 2 7r 



(90) 



(91) 



It is interesting to note that thanks to the renormaliza- 
tion group equation pOjl. we can directly compute Ag T 
at the scale fio and only with the help of Ea. ((81|l . we have 
access at the scale fi. 



A. One loop energy Regularization and 
Renormalization for the trace 



As for the Eq. l(72|) , the energy density associated with 

^trace 



A' 



2 1 



3 64tt 2 



u} 2 \/uj 2 — ^„ du 



2 to 4 (r) 

3 1024tt 2 



m 2 (r) 



In ( -4 ! +21..2- - =p*7/ e ( M ,r) 

(98) 

As for the quantity in Eq. I|79|l , we have a dependence on 
the arbitrary mass scale \x. The use of the renormalization 
group equation gives 

1 d t 



8ttG 

Solving Ea. H99|) . we find a running A race provided that 
the scale fj, is changing 



A trace 
Ji 



(/i,r)=A ™ ce ( Mo ,r) + 



Substituting Ea.l(TUD)l into Eq.(0U we find 

A race (// ,r) 
8ttG 



2 1 



-foo 



3 64tt 2 j 
where we have defined 



w 2 \/w 2 -m 2 (r)duj, (92) 



to 2 (r) 



MG 



(93) 



The use of the zeta function regularization leads to 



2 m 4 (r) 

3 1024tt 2 



hi 



m 2 (r)\ 1 



(101) 



If we adopt the same procedure of finding the minimum 
for A ™ ce (/J,Q,r) we discover that the only consistent so- 
lution is for M = 0. This leads to the conclusion that 
the trace part of the perturbation does not contribute to 
the cosmological constant. 



p trace ^ 



2 1 



3 64tt 2 



+ 00 

■ 2 I dtu- 
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(94) 



In order to renormalizc the divergent trace contribution 
to ZPE, we observe that from Eq. (|72fl . we get 



A trace = ^Qptrace 



(95) 



To handle with the divergent energy density we extract 
the divergent part of A <roce , in the limit e — > and we 
set 



A 



trace, div 



2 G 



3 128tt£ 



m 4 (r) . 



(96) 



In this paper, we have considered how to extract 
information on the cosmological constant using the 
Wheeler-De Witt equation. In particular, by means 
of a variational approach and a orthogonal decompo- 
sition of the modes, we have studied the contribution 
of the transverse-traceless tensors and the trace in a 
Schwarzschild background. The use of the zeta function 
and a renormalization group equation have led to Eq. I|91|) 
and recalling Ea. l|39|> . we have obtained 



M, f) 



A trace 
^0 



(M,i 



(102) 
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If we choose to fix the renormalization point fiQ = m p , 
we obtain approximately Aq (M, rj ~ 10 37 GeV 2 which, 
in terms of energy density is in agreement with the es- 
timate of Eq.lQ. Once fixed the scale po, we can see 
what happens at the cosmological constant at the scale 
p, by means of Eq.JHIJ. What we see is that the cosmo- 
logical constant is vanishing at the sub-planckian scale 
p = m p cxp (— t), but unfortunately is a scale which is 
very far from the nowadays observations. However, the 
analysis is not complete. Indeed, we have studied the 
spectrum in a W.K.B. approximation with the following 
condition kf(r,l,uji) > 0, i = 1,2. Thus to complete 
the analysis, we need to consider the possible existence 
of nonconformal unstable modes, like the ones discov- 
ered in Refs.0. If such an instability appears, this does 
not mean that we have to reject the solution. In fact 
in Ref . | 111] - we have shown how to cure such a problem. 
In that context, a model of "space-time foam" has been 
introduced in a large N wormhole approach reproduc- 
ing a correct decreasing of the cosmological constant and 
simultaneously a stabilization of the system under exam- 
ination. Unfortunately in that approach a renormaliza- 
tion scheme was missing and a W.K.B. approximation on 
the wave function has been used to recover a Schrodinger- 
like equation. The possible next step is to repeat the 
scheme we have adopted here in a large N context, to re- 
cover the correct vanishing behavior of the cosmological 
constant. 



APPENDIX A: THE ZETA FUNCTION 
REGULARIZATION 

In this appendix, we report details on computation 
leading to expression (J73J. Wc begin with the following 
integral 

,2 E r+°° 



I + = p^C°°du> ^ r 



(Al) 



with m 2 (r) > 0. 



1. J+ computation 



If we define t = u>/ \Jrv? (r), the integral 1+ in Eq. (|A1|) 
becomes 



p(e) = p 2£ m l -* e {r) / dt 



t 2 



o {t 2 + If 



1 2e 4-2=/ N r (l) r ( e - 2 ) 

—p ui Ze (r) —^4 — 

2 r( e -i) 



4 (n \m 2 {r)) r (e-f)' 



(A2) 



where we have used the following identities involving the 
beta function 



B(x,y) = 2 dt- 
Jo ( 



/o (t 2 + l) x+y 

related to the gamma function by means of 

T(x)T(y) 



Rex > 0, Re y > 
(A3) 



B(x,y) = 



T(x + y) 



(A4) 



(A5) 



Taking into account the following relations for the T- 
function 

r ( e - 2) = r(1+e) 

1 v fc A ) - e (e-l)(e-2) 

and the expansion for small e 
r (1 + e) = 1 - 7 e + (e 2 ) 

r (^ + 5)= r(i)- e r(i)( 7 + 2in2) + o( £ 2 ) 

x £ = 1 + e\nx + O (e 2 ) 



where 7 is the Euler's constant, we find 
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(A6) 
(A7) 



2. 7_ computation 



If we define t — u/ \J m 2 (r) , the integral 2_ in Eq. (|A1|) 
becomes 



p (e) = p 2e m 4 ~ 2£ (r) / dt 



t 2 



2 T\ £ ~2 



(t 2 -l) 



l ~^m^{r)B[e-2,l-e 



r(f- g )r( g -2) 



' rn*(r)( z £ : Xr(*-:)T { :-2). (As; 



4v^ 



m 2 (rj 



where we have used the following identity involving the 
beta function 



i B (i-v-<±,v)=j 1 + °°dtt»-Htp-iy 



p > 0,Ref > 0,Rep, < p — pRev 



(A9) 
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and the reflection formula 

r 0) r (i - z) = -zT (-z) r (z) 



we find 



(A10) 



From the first of Eas. (|A5(l and from the expansion for 
small e 

r (f- £ ) = r (|) (l-£(-7-21n2 + 2)) + 0( £ 2 ) 



m 4 (r) 
16~ 



1 + e In x + O (e 2 



(All) 
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m 2 (r) 



2 In 2 



(A12) 
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